Thermal fluctuations and anomalous elasticity of homogeneous nematic elastomers 
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We present a unified formulation of a rotationally invariant nonlinear elasticity for a variety of 
spontaneously anisotropic phases, and use it to study thermal fluctuations in nematic elastomers 
and spontaneously anisotropic gels. We find that in a thermodynamic limit homogeneous nematic 
elastomers are universally incompressible, are characterized by a universal ratio of shear moduli, and 
exhibit an anomalous elasticity controlled by a nontrivial low temperature fixed point perturbative 
in D — 3 — e dimensions. In three dimensions, we make predictions that are asymptotically exact. 



Qualitatively important effects of fluctuations are 
generically confined to the vicinity of isolated critical 
points, where a system is tuned to be "soft", charac- 
terized by low energy excitations. As long as the or- 
dered state is stable, fluctuations about it are typically 
described by a harmonic theory, controlled by a Gaus- 
sian fixed point. However, there exists a novel class 
of systems, that includes smectic[[j], || and columnar 
liquid crystals H. vortex lattices in putative magnetic 
superconductors [B , polymerized membranes and ne- 
matic elastomers^, 0, |§, whose ordered states are a 
striking exception to this rule. A unifying feature of these 
phases is their underlying, spontaneously broken rota- 
tional invariance, that strictly enforces a particular "soft- 
ness" of the corresponding Goldstone mode Hamiltonian. 
As a consequence, the usually small nonlinear elastic 
terms are in fact comparable to harmonic ones, and there- 
fore must be taken into account. Similarly to their ef- 
fects near continuous phase transitions, but extending 
throughout an ordered phase, fluctuations drive nonlin- 
earities to qualitatively modify such soft states. The re- 
sulting strongly interacting ordered states at long scales 
exhibit rich phenomenology such as a universal nonlocal 
elasticity, a strictly nonlinear response to an arbitrarily 
weak perturbation and a universal ratio of wavevector- 
dependent singular elastic moduli, all controlled by non- 
trivial infrared stable fixed points illustrated in Fig.[l]. 
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FIG. 1: Renormalization group flow for anomalously elas- 
tic solids, with T c the transition temperature to the ordered 
state and T a nontrivial infrared stable fixed point controlling 
properties of the strongly interacting ordered, critical state. 



clastic Goldstone mode corresponding to a spontaneous 
orientational ordering and an accompanying elastic dis- 
tortion of the network H pi [l2|]. Despite of considerable 
progress in understanding the properties of these liquid 
crystalline rubbers[Q, it was only recently appreciated, 
that (because of their harmonic softness) nematic elas- 
tomers and gels belong to the aforementioned class of 
anomalously elastic materials, and are the first example 
of three-dimensional solids with such novel fluctuation- 
driven properties. 

In this Letter we present a unified formalism for a sys- 
tematic derivation of the rotationally-invariant nonlin- 
ear elastic energy for variety of soft elastic media. This 
is a nontrivial ingredient that is a prerequisite to treat- 
ing effects of fluctuations and external perturbations in 
such soft materials. We utilize this formalism to study 
long-scale elastic properties of homogeneous fluctuating 
nematic elastomers, and demonstrate that they indeed 
exhibit anomalous elasticity phenomenology discussed 
above || . 

A configuration of an elastic media is described by a 
mapping x t— > i?(x) from the D-dimensional reference 
space x, which labels mass points of the media, to the 
d-dimensional target space R. Our goal is to construct 
a nonlinear elastic energy functional TL, describing long- 
scale phonon deformations u(x) = i?(x) — i?°(x) about 
the spontaneously ordered state i?°(x), with TL exhibiting 
all rotational and translational symmetries of the refer- 
ence and target spaces |^]. This is nontrivial to imple- 
ment at a nonlinear level as some of the symmetries are 
spontaneously broken by the ground state, and are there- 
fore not manifest (hidden) in TL. 

A systematic way of constructing such TL is by expand- 
ing it in in terms of scalar nonlinear operators S n 



Nematic elastomers and gels are weakly crosslinked 
polymer-liquid crystal composite materials, which, when 
nematically ordered exhibit fascinating elastic and 
electro-optic properties, and are therefore of considerable 
technological importance. The most unusual of these is 
the divergent strain response to a stress applied trans- 
versely to the spontaneous nematic direction Q, an 
elastic softness that is a consequence of the novel nemato- 



n[g] = TL[g 



S n = Tr(g»-£") 
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with 



9ij — diR ■ OjU ana 
tensors corresponding to -R(x) and -R°(x) configurations 
and summation convention over repeated indices is used 
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throughout. Eliminating the metric tensor 
9*3 = 9ij +2-^ zr u a p- 



dxi 



in favor of the Lagrangian strain tensor 
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measured in terms of the ground state coordinate sys- 
tem (rather than x{) leads to the desired nonlinear 
elastic energy Ti[u a p\. The necessity of this procedure is 
that it ensures that the expansion in u a p about the spon- 
taneously symmetry-broken ground state is rotationally 
invariant in both the target and reference spaces. 

We have used this general formalism, with differences 
encoded in £>, d and the nature of the broken ground 
state, R° (x) , to derive nonlinear elastic Hamiltonians for 
a variety of elastic media [fl4|, some of which have previ- 
ously appeared in the literature. In this Letter wc apply 
it to study fluctuations in the most nontrivial and hereto- 
fore unexplored system, the spontaneously-uniaxial ne- 
matic elastomers with d = Z?j~ 

As discussed in detail in Refs. 
erties of nematic elastomers can be captured by a purely 
elastic description in terms an elastic strain tensor, with 
a nematic order parameter Q a p — S(n a fip — 5 a p/3) (S 
and n, respectively, the magnitude and the unit director 
for nematic order) integrated out. The isotropic-nematic 
transition is then characterized by a spontaneous uniaxial 
distortion u^p oc Q a p, corresponding to a ground-state 
conformation and a metric tensor given by 




most of the prop- 



i?°(x) = (±X±+(zZZ, 
9ap = Qa°aP> 



(0) 
(7) 



with x^ a D — 1-dimcnsional vector transverse to the 
nematic director n = z. 

In the physical case of D = d = 3 there are three non- 
linear rotationally invariant operators S n , which, using 
Eqs.| and |, can be readily shown to be given by 

Si = 2(lu aa , (8) 
<S 2 = 4C„u QQ + 4:(l(pu a pup a , (9) 
S3 = 6( a u aa + 12C,aCpU a pup a + SClCpC^u a pup 7 u ja , (10) 

Using these expressions inside the expansion for TL 1 
Eq.^, taken to quadratic order in S^s, we obtain 
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where Roman indices (i,j,k,...) take values in the _L- 
space only and 
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l^zi 



2 ai Cl ± +4a 2 CU + 6a 3 (i ± , (12) 
16a 2 C z 2 C! + 24a 3 (C! + CzKzCl, 



(CX-cV), 



(13) 



with the detailed form of other elastic constants not im- 
portant here|l^|. 

In the absence of fluctuations, our expansion about 
the ground state g a , characterized by ensures that 
o-z = a_L = in equilibrium. Given that in the nematic 
ground state, Cz 7^ C-L> this then leads to a strict van- 
ishing of the shear modulus jx Z i = 0, as anticipated by 
Golubovic and Lubenskyjf^, |8| and is implicit in the neo- 
classical theory of nematic elastomers^, |llj. It corre- 
sponds to a vanishing energy cost of a shear distortion u Z i 
displayed in Fig.|^a, that is one of many fascinating prop- 
erties of nematic elastomers and is a manifestation of the 
"soft" Goldstone mode elasticity, discussed in the intro- 
duction. At finite temperature, fluctuations renormalize 
all elastic constants and in particular generate nonzero 
a z , a±. However, Ward identities associated with the 
underlying rotational invariance of TL ensure that such 
linear in strain terms can always be eliminated by shift- 
ing to a true, thermally renormalized spontaneous defor- 
mation g° Q. 

Keeping only the most relevant terms in 7i, Eq.[l2], and 
using nontrivial relations between elastic coefficients p], 
enforced by the rotational invariance, we arrive at the 
effective nonlinear elastic Hamiltonian for a uniaxial ne- 
matic elastomer: 



Hne = -^(w zz +wu) 2 +C(w zz +wu){w z 
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where 



w. 



= d z u z + -(V_lu z ) 2 , 

= IjidiUj + djm) - ^diu z dju z 
w k k , 



D - 1 



(15) 



are rotationally-invariant smectic-like and columnar-like 
nonlinear strain tensors relative to the uniaxial state, 



B = ^(B z +B ± + X zl ), 



fJ-L 



C 



(B z + B± — X zi ), 
{B z - B ± ), 



(16) 
(17) 
(18) 



and B± = A + 2/i/(Z? — 1) is the _L-plane bulk modulus. 
It is easy to see that B is the overall bulk modulus, \il 



3 



and p are the longitudinal and transverse shear moduli, 
respectively corresponding to deformations displayed in 
Fig. U b,c, and C couples bulk mode with longitudinal 
shear. Symmetry requires that transverse shear mode 
remains decoupled from other modes. Because of the 
vanishing jjL Z i, we have also added a curvature K elastic 
term for u z to ensure stability at long scales. 
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FIG. 2: (a) simple, (b) longitudinal and (c) transverse shear; 
in (c), nematic order is out of the paper. Simple shear trans- 
forms the system into an equivalent nematic ground state, 
and thus costs no energy. 

Using Hne, we now proceed to study effects of thermal 
fluctuations on nearly homogeneous nematic elastomers. 
Standard analysis leads to the Fourier transform of har- 
monic correlators G^Jx) = (u Q (x)w | g(0))o: 



GL(q) 

G°,(q) 
G°(q) 



1 



B z {l-p) q l + Kq\ 
( pB z 



\2p + \J ql{B z (l- p)ql+Kqiy 

B z q 2 + Kq\ q iqj 

{B z (\-p)ql + K q \) (2 M + A)gl 

PQ± 



(19) 

(20) 
(21) 



where p = (B — p L ) 2 /{B z {2p + A)). 

To assess the role of nonlinear elastic terms in the 
presence of thermal fluctuations we compute perturba- 
tive corrections to elastic constants in TIn e and find that 
they are dominated by u z fluctuations. Given the struc- 
ture of G° zz correlator, as in smectic liquid crystals this 
perturbation theory diverges for D < 3 on length scales 
longer than ^ L = A- 1 g~ 1/{3 - D) (= A~ 1 e 32 / ( - 59g± \ for 
D = 3), with g± defined in Eq.^4[ below, and A an ultra- 
violet cutoff associated with the elastomer mesh size. 

To understand the elastomer properties beyond £,nl, 
we perform a momentum-shell renormalization group cal- 
culation, perturbatively in e = 3 — D. Integrating out 
short scale phonon fluctuations u a and rescaling the spa- 
tial coordinates and fields, so as to keep A fixed, we arrive 
at the following flow equations for the elastic constants 
at scale A~ 1 e e : 
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(22) 



;9±, 



(23) 



Vk 



9L9± (9 + Upx - 22^p 2 ) - 4g L 2 (-1 + p 2 2 
8 (g± Pi+9l (1 + Pi - 2 y/plp 2 )) 



The physics is controlled by the flows of two dimen- 
sionless coupling, gh{£) & n d g±(£), and two dimensionless 
ratios p\{€) and p 2 {t): 



9± 



/' 



(B-2C + p + p L ) 



4nKi V Bp + 2C (-2C + p) + (4B+/i) p L 

C 



9l = — 9x, Pi = -5-, P2= m (24) 

p B \JBp h 

Flow equations for p\{€) and p 2 (£) are given by: 

^ = (25) 

£ = 4«»<i-A, (26) 

where mechanical stability requires \p 2 \ < 1- With the 
exception of an unstable fixed point at g* L — and an un- 
stable fixed line at \p 2 \ = 1 (parameterized by pi)fil|], we 
expect (and verify a posteriori) that the infrared stable 
fixed point is characterized by g* L 7^ 0. Equations [2^j2t^ 
then imply that at such fixed point, p\ = p 2 = 0. The 
flow equations for g^ (£) and g± (£) , the full form of which 
is too involved to be reproduced hereQ, then simplify 
considerably: 



dgL 
dl 

dg± 
d£ 



£9l 



9L (40ff L 2 + 68g £g _L + 13gj 

8(4g L +g ± ) 
ff-L (4.9l 2 + 32. 9L g ± + 7g ± 2 
4(4g L +g±) 



-, (27) 



(28) 



In addition to reproducing the 3 — e-dimensional thermal 
smectic fixed point [Q, with g\ = 0, g* L = 4e/5, we find an 
infrared stable fixed point at g* ± — 32e/59,g2 = 4e/59, 
with 

TIB = 0, r^ = 38e/59, (29) 
Vc = Vl = Vl =ri±= 4e/59, (30) 

characterizing long-scale anomalous elasticity of homoge- 
neous nematic elastomers Iful. 
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For the physically most interesting case of D = 3, 
and g±(£) are marginally irrelevant, allowing an 
asymptotically exact computation of anomalous elastic- 
ity. Standard asymptotic analysis of Eqs. p5f|28| with 
e = 0, shows that for large I coupling constants flow 
to zero algebraically: 



P2(i) 

9l{Z) 



£-4/59 ] 
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(31) 
(32) 

(33) 



Matching calculations, together with Eqs.|22| then pre- 
dict that on long scales all effective elastic constants 
except B are anomalous, depending logarithmically on 
wavevector q of the deformation 



B(q) <* Bo, (34) 
.C(q) oc (log|q ± |)- 4 / 59 , (35) 
K(q) oc (log|qJ) 38 / 59 . (36) 



Hence at long scales, the bulk modulus is much larger 
than shear moduli, predicting that independent of micro- 
scopic details nematic elastomer is effectively incompress- 
ible in the thermodynamic limit. Logarithmically diver- 
gent strain-stress response immediately follows. Equa- 
tions |24|, predict a universal ratio of the renormalized 
shear moduli 



(37) 



Although detecting logarithmic length scale dependence 
predicted here is likely to be difficult , the incompressibil- 
ity of nematic elastomer and above universal ratio should 
be readily observable. 

Based on our predictions in 3 — e dimensions, we would 
expect a much stronger, power-law anomalous elastic- 
ity in two-dimensional (2D) nematic elastomers, with 
anomalous exponents given by Eq.|9|j30| withe = 1. How- 
ever, a moment of reflection shows that because in 2D the 
subspace perpendicular to the nematic direction is one- 
dimensional, there is no transverse shear mode. We must 
therefore set [i = 0, or deduce the properties 2D elas- 
tomers from a D — 1-axial analytical continuation, as we 
have recently done for nematic elastomer membranes Jl5|. 
A more serious difficulty is that in our analysis we have 
ignored the subdominant columnar-like elastic nonlincar- 
ities, that, although irrelevant near D — 3, are relevant 
for D < 5/2 and in two dimensions are in fact identical to 
the smectic-like nonlinearities that we have studied here. 
Simultaneously treating both nonlinearities remains an 
interesting and challenging problem. 

Throughout this work, we have also ignored elastomer 
heterogeneities. However, real elastomers are only sta- 
tistically homogeneous and isotropic, and therefore re- 
sults presented here are only valid up to a long but finite 



disorder-dependent length scale £^ L . On longer scales 
network heterogeneity dominates over thermal fluctu- 
ations and results presented here crossover to a dis- 
tinct phenomenology, controlled by a nontrivial zero- 
temperature fixed point ||, ||, |^, |). 

To summarize, we have presented a unified formulation 
of a rotationally invariant nonlinear elasticity for a vari- 
ety of elastic media with spontaneously broken spatial 
symmetry. We applied this formalism to nematic elas- 
tomers and studied their thermal fluctuations. We found 
that they exhibit anomalous elasticity, which in three 
dimensions is characterized by a logarithmic wavevec- 
tor dependence of elastic moduli, universal ratio of shear 
moduli, and a long-scale incompressibility. 
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